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Abstract. We propose a simple model with two infective classes in order to model the 
cholera epidemic in Haiti. We include the impact of environmental events on the epidemic 
in the Artibonite and Ouest regions by introducing terms in the force of infection that vary 
with environmental conditions. We fit the model on weekly data from the beginning of the 
epidemic until March 2012. We then used this model to obtain epidemic projections from 
April 2012 through February 2013, and then modified these projections incorporating the 
vaccination programs that were recently undertaken in the Ouest and Artibonite regions to 
compare with actual cases. 

Using real-time daily rainfall we found lag times between precipitation events and new 
cases that vary seasonably, ranging from 5 to 11 weeks in Artibonite, and 5.7 to 8.5 in Ouest. 
In addition, it appears that, in the Ouest region, tidal influences play a significant role in the 
dynamics of the disease. 

Intervention efforts of all types have reduced case numbers in both regions, however, per- 
sistent outbreaks continue. In Ouest where the population at risk seem particularly besieged 
and the overall population is larger, vaccination efforts seem to be taking hold more slowly 
than in Artibonite where a smaller core population was vaccinated. With the implemented 
vaccination program the model projects that by the end of February 2013, the number of 
cases in Artibonite would be decreased by about 7 thousand, and in Ouest by about 3 thou- 
sand. 

Keywords: cholera, Haiti, bioclimate, epidemic model, vaccination. 

I. Introduction 

On January 12, 2010, a 7.0 magnitude earthquake struck near Haiti's capital, Port-au- 
PrinceP The poorest nation in the Western Hemisphere, the earthquake shattered Haiti's 
already weak infrastructure.^ Thousands of Haitians were killed and even more were forced 
to flee to resettlement campsP 

In October 2010, the first case of cholera was reported in Haiti. A later UN investigation 
revealed the specific strain of V. cholerae came from South AsiaP The UN investigation and 
epidemiological literature suggest that the epidemic began outside of a UN peacekeeper camp 
near Mirebalais in the Centre department, along the Artibonite River EH3 As V. cholerae is a 
waterborne pathogen, the Artibonite River is the ostensible route through which the disease 
spread throughout Haiti's ten administrative regions, called departments™ 

Anecdotal news reports describe the dismal situation for thousands of Haitians that still 
remain displaced months after the earthquake. Sewage of millions of people flow through 
open ditches. Human waste from septic pits and latrines is dumped into the canals, and after 
it rains, ends up in the sea. Those living close to the water use over-the-sea toilets, and next 
to these outhouses, fishing boats unload and sell the fish from plastic buckets.. P 

Haiti's two most populous regions, Ouest and Artibonite, were also the two regions hardest 
hit by the epidemic. Cases in Ouest and Artibonite account for 60% of the total burden of 
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cholera in HaitiP For this reason, we chose to focus our analysis on the Ouest and Artibonite 
regions. By April 7, 2012, cholera had affected 5.7% of the total population^] in Ouest and 
6.9% of the population in ArtiboniteP 

We reviewed the literature dealing with cholera, modeling, and climatic conditions (rainfall, 
precipitation, & tides) in Haiti at the beginning of this study and again at the end. We did 
multi-database searches in the Biology, Medicine, and Health fields and found four modeling 
papers^H^ of which the first three used variations of a basic system of differential equations 
proposed in 2001 by CodegoJm At that time, none of these models took environmental 
conditions directly into account. Two other papersSHGH have since appeared that do take 
precipitation into account in cholera. The firslP^ is a spatiotemporal Markov chain model 
using seasonal rainfall which was outside our purview. The second^ deals specifically with 
Haiti and was done by the same group that produced one of the earlier papersP Here they 
looked at the reliability of these earlier studies, and they found that although these models 
do well in capturing the early dynamics of the epidemic, they fail to track latter recurrences 
forced by seasonal patterns.^ As a follow up, Rinaldo et al.^ add a precipitation forcing 
function to their original model along with other modifications such as the river network and 
population mobility. These modifications produce a better fit to the observed pattern of case 
but predict recurring large outbreaks tracking the seasonal precipitation patterns. 

The three models^' proposing a variation on the SIWR model proposed by Codego in 2001^ 
explain cholera transmission through susceptibles' contact with a water reservoir, rather than 
susceptibles' contact with infectious individuals. The models proposed by Tuite, et alP and 
Bertuzzo, et alP both incorporated a "gravity" term to study the interaction among regions. 
The model proposed by Andrews and BastP accounted for a bacterial "hyperinfectivity" stage, 
following research by Hartley, et al. in 2006 showing that V. cholerae initially has a higher 
infectivity before it decays to a lower infective rate in the aquatic reservoir. The fourth model 
by Chaoa, et alP^ was an agent based approach designed to look at various vaccination 
strategies. 

In addition, all these models assessed the impact of potential intervention strategies, in- 
cluding vaccination. Bertuzzo found that a vaccination campaign aiming to vaccinate 150,000 
people after January 1, 2011 would have little effect, in part because of the late timing and 
in part because of the large proportion of asymptomatic individuals who would need to get 
vaccinatedP Both the models proposed by Tuite, et alP and Andrews and BastP suggest 
that vaccination campaigns would have modest effect. In March 2012, Partners in Health 
began vaccinating 100,000 individuals with Shanchol, a two-dose cholera vaccine.^ The size 
of the campaign is limited by the size of the global stockpile of Shanchol.^ The vaccination 
campaign is targeted at 50,000 individuals living in the slums of Port-au-Prince, where pop- 
ulation density is thought to increase the rate of cholera exposure, and at 50,000 individuals 
living in the Artibonite River valley, where the epidemic beganP^ Chao et alP^ showed that 



*Note: population data for Haiti is from 200!^ , one year before the earthquake. 
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a targeted vaccination strategy would have the best results for this limited supply of vaccine, 
and by vaccinating 30% of the population the cases could be reduced by as much as 55%. 

In 2001, Codego proposed introducing an oscillating term to model seasonal variability.^ 
However, none of the four Haiti-specific models accounted for seasonality. Haiti experienced 
flooding in June 2011, October 2011, and March 2012.^1 As cholera has reached an en- 
demic state in Haiti, an analysis of cholera's seasonality as it relates to Haiti's rainy season 
is pertinent. Moreover, mathematical models should incorporate seasonality in order to more 
accurately predict the course of the epidemic and to simulate the effects of potential inter- 
ventions. In April 2012 Rinaldo et alP^ reexamined the above four models (including their 
ownP) and concluded that, among other factors, seasonal rainfall patterns were necessary to 
account for resurgences in the epidemic. They use long-term monthly averages to augment 
the bacterial growth term of contaminated water-bodies. 

Other papers dealing with bio-climate were a study of Cholera in Zanzibar, East Africa 
that demonstrated 8 weeks fixed delay^ between rainfall and cholera outbreaks, and another 
paper in Bangladesh^ that have a somewhat shorter delay (4 weeks). Both these studies use 
a statistical approach with seasonal data. Both also made note of the potential influence of 
ocean environmental factors, and the Reyburn et al. paper included sea surface height and sea 
surface temperature in their analysis but failed to find any significant relationship. 16 In a third 
paper Koelle et al. (2005) 13 model very long time periods, more than a year, in Bangladesh. 
This model also uses seasonal precipitation and models changes in the susceptible fraction of 
the population due to demographics and loss of immunity. 

In this paper we use detailed and current rainfall, temperature, and tidal records to model 
cholera in the Artibonite and Ouest regions. We forego a bacterial compartment in favor of 
an a posteriori approach using climatic data directly to estimate infection rates. This has the 
advantage of more tractable temporal estimates without over parameterizing and including 
compartments that are essentially unmeasurable. This paper is also the first, that we know 
of, that uses tidal variation in a model of cholera dynamics. 

These long term trends and environmental influences established the pattern of response of 
the epidemic in Artibonite and Ouest. Thus parameters were chosen and model calibration 
set prior to a vaccination program being implemented. We then used the model to evaluate 
the performance of the vaccination program against the backdrop of an alternative history 
without vaccination. 



2. Material and Methods 

For Ouest and Artibonite, we investigated the correlations between reported cholera cases 
and rainfall, temperature, and in the case of Ouest, tidal range. We wanted to determine a) 
if such correlations exist, b) the time delay between environmental conditions and recorded 
cholera outbreaks, and c) if effectiveness of a recent vaccination program could be assessed 
by use of this model. 
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2.1. Model. We take a combined mechanistic and phenomenological approach to our model. 
The mechanistic part is a standard SIR type model where individuals move from Susceptible 
(S) to Infected to Recovered (R) classes^ An infectious individual may either be symptomatic 
(I) or asymptomatic (A)W^ The probability, p, of asymptomatic infection is 0.79J22E1 Both 
symptomatic and asymptomatic individuals move to the recovered group, R, at a rate 7. 
Symptomatic individuals die from cholera at a rate p (dead are denoted by D). 

The phenomenological part comes from estimating the force of infection, /3(t), by fitting 
the number of cases predicted by the model (both incidence and cumulative cases) to the 
data. This has the benefit and drawback of eliminating the connection, in the model, between 
the infected population and the infecting water bodies, food, and so on. A benefit since 
we don't need to model and estimate the extra compartments and half dozen or so other 
parameters necessary to incorporate the mechanism of environmental bacterial populations 
and their containing water bodies. And a drawback because an important mechanism is no 
longer included in the model. We are obviously not suggesting that such feedbacks don't 
exist, only that they are beyond the scope of this study. All other parameters are discussed 
below and values are given in Table [T] and Table [2] 

This paper evaluates Artibonite and Ouest separately, in order to capture the varying 
dynamics in each region. We use the following system of ordinary differential equations, and 
the components to be fitted: N = S + A + 1 + R + D, where N is the initial population of the 
region. (For the purposes, and time scale, of this study we chose to ignore the demographics 
of the population.) 

System of equations Variables to be fitted 

g = -/?(*) S ) 

I = pP(t) S- 1 A I f new cases = (1 - p) Jl At 0(0 S(0 d£ \ 

f = (1 - p)/3(t) S - (7 + y)I [ I total cases = (1 - p) f< (3(0 S(£) d£ j 
l = 7(A + /) J 

2.2. Environmental Effects. We assume that significant rain events will cause flooding, 
which will increase direct contact the bacteria via contaminated water, or indirectly via soils, 
vegetation, and disease bearing insects, etc. that have been contaminated or consumed cholera 
bacteria. We also assume that temperature plays a dual role by increasing the infection rate 
and decreasing the lag time between rainfall events and disease outbreaks. 

Additionally, in the Ouest region there are commercial areas, tent-cities and slums that have 
raw waste directly discharging to Port-au-Prince Bay, or to the bay via rivers (e.g. Froide, 
Momance, and Grise), and numerous open sewage canals. We hypothesis that there may be 
additional disease generated when large tidal ranges stir up contaminated sediments or cause 

'Conceptually the underlying model is a variation on the SIWR model, which assumes that cholera is 
spread through susceptibles' contact with contaminated water, food or fomites. This model uses the amount 
of water consumed as a proxy for all possible modes of transmission, and the concentration of bacteria in the 
water consumed modifies the infection rate by a dose-response expression (see^ with a base mode^^) 
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blooms of plankton in these coastal waters. Again the chain of infection may be complex and 
include direct contact with water, insects, or consumption of contaminated seafood, etc. 
Therefore, we modify the force of infection term, /3, such that 

(1) p(t) = u(t) \a p H{t) P(t - Tp , 6 P ) + a m M(t - r m , d m )\ , 

where P is a moving average of the amount of daily rainfall, and M is a moving average of 
the maximum semi-diurnal tidal range, if is a heat index based on mean air temperature, r p , 
is the lag time for precipitation (which is itself affected by temperature), r m is the lag time 
for tides, 9 P , and 9 m are the averaging periods, and a p , and a m are proportionality constants. 

We also included an expression for improvement of conditions over time, either through 
reducing the effective susceptibles, or reduction in the infection rate, or both. These fac- 
tors would be due to increased access to clean water, increased personal hygiene, decreased 
contamination of the environment, vaccination, and rapid treatment of new cases, or any 
other means of removing risk. We assume the force of infection is reduced by an exponential 
decreasing factor u(t) = e~ rt + (1 — e~ rt )uo, where r is the rate of improvement in conditions 
and 1 — Uo is the maximum amount of improvement possible in the short term (uq can be con- 
sidered the persistent fraction of risk remaining in a region). This approach is tantamount to 
having the susceptible, or at-risk population reduced directly by public health improvement^] 

2.3. Data. The source of the epidemic data was the Haitian Ministry of Public Health and 
Population and compiled by the Pan American Health Organization^^ The data sets avail- 
able, which we used for this study, consist of cumulative cholera cases, and new cholera cases. 
New cholera cases are calculated based on the difference between the latest report and the 
previous one. Cholera case definition also includes suspected cholera cases and deaths in 
addition to confirmed cases and deaths. As such, data posted on the web site are periodically 
updated with minor corrections. Cases are reported on weekly basis with the reporting week 
beginning on Sunday. 

Reported hospitalized cases and hospitalized deaths are probably more accurate, but for 
the purposes of this study, less useful, since we want to track the progress of the epidemic, 
and increases in access to treatment biases the data. New case estimates for the first few 
weeks of the epidemic are unavailable. Total cases (TC) for the first 4 weeks was obtained by 
regressing, for the following 8 weeks (14-Nov-10 through 2- Jan- 11), the reported total cases 
against hospitalized cases (HO). For the Ouest region TCo = 2.2HCo (R 2 = 0.99) , and 
for the Artibonite region TC A = 2.59HC A + 497 (R 2 = 1.00) . 

Rainfall, temperature and tide data were compared to the pattern of new cholera cases 
reported each week. The rainfall data comes from NASAp^ using centrally located points 
in each region as our datum point, given by (latitude, longitude). For Artibonite our da- 
tum point is (19.125,-72.625) and for Ouest, it is (18.625, -72 • 375). 22 NASA updated 

tFor example /?(t) = a p H(t) P(t - r p , 6 p ) + a m M(t - r m ,6 m ) , § = -f3(t) S - r (l - ^) S, and f = 
— r (u — Uo) , etc. Where the rate of reduction of risk of the susceptible population starts at r (1 — u ) and 
decreases to when u = Uq. 
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its processing of the data during the course of this project from TRMM_3B42_daily.006 to 
TRMM_3B42_daily.007. Although our initial investigation used the version 6 material, all 
figures and parameter estimates given in this paper are using version 7. Precipitation esti- 
mates are a combination of remote sensing and ground verified information reported with a 
spatial resolution of 0.25 x 0.25 degrees and a temporal resolution of 1 day, further details are 
available on the web sitepZl 

Temperature data is mean daily air temperature at Port-au-Prince, and is reported by 
the Weather Underground.^ Although we had daily numbers for Temperature we chose to 
extract the annual cycle and use that in simulations. This has the advantage of allowing us 
to extrapolate temperature for model projections. 

Tide data is for Port-au-Prince Bay (Stationld: TEC4709). 24 These numbers are predic- 
tions from NOAA's tide model for this location and are not direct measurements. NOAA's 
web site allows you to download tide numbers for any date starting from 2010 and extending 
through 2014. 

The initial modeling was done by comparing data sets in the frequency domain for new cases 
and rainfall in order to find any suggestions of matching periodicity and/or time-lags. Code 
was then written in Berkeley Madonna to simulate the dynamical system, and try various 
manipulations of the environmental data in order to match predictions to data^J Rainfall data 
was available only to week of 30-Sept-2012 (week 102) because of the lag between the timing 
of rainfall events and what had been processed at the time we accessed itP^l 

2.4. Model calibration. All modeling and calibration was done in Berkeley Madonna. 
When possible we used reasonable ranges for each parameter as established by previously 
published research (see Table |2|. Since our model contains several fitted parameters, we at- 
tempted to minimize the plausible range for each parameter in order to obtain a more accurate 
fit. Parameter values not obtained from the literature were adjusted in an iterative fashion to 
narrow their initial ranges, fitting by eye predicted number of new cases to the reported num- 
ber of new cases. The Berkeley Madonna curve fitting algorithm was then used to minimize 
the root mean square difference between model predictions of number of cumulative cases and 
reported cases for the restricted parameter space. Confidence and prediction intervals were 
calculated using the delta method.^ 

2.5. Parameters from literature. Table [T] displays parameter values for each region ob- 
tained from published or online sources. 

2.6. Environmental components of the force of infection. 



f When we discuss the output of our model, will use the term "prediction" to indicate the numbers generated 
during the calibration phase of the model and "projection" to indicate the extrapolation of the model past 
cases used for the calibration. 
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Parameter 


Artibonite 


Ouest 


Units or calculation 


References 


n 

r 


0.79 


0.79 


fraction becoming asymptomatic 


|20||26| 


7 


1.4 


1.4 


fraction recovered per week 




N 


1,571,020 


3, 664, 620 




® 


So 


1,534,338 


3,663,699 


N - {A + Jo + Ro + D ) 




h 


7,653 


193 


regression on early hospitalized cases 


ED 


A 


28, 790 


726 


rVo 




Ro 












Do 


239 


2 






to 


17-Oct.-2010 


17-Oct.-2010 


290 th day of the year 


m 


Uq 


0.05 


0.05 


persistant fraction of pop. at risk 


ED 



Table 1. Parameter values obtained from literature. 



2.6.1. Artibonite. Bioclimatic influences on infection rate are a product of precipitation rate 
and an heat index. Daily precipitation p (tf^ is averaged over an interval 6 P so the running 
average precipitation rate is 

p(t-T p ,e p ) = T Y,p(t-T P -j), 

where t p is the delay in precipitation's affect on infection rate (see below). For the delay 
we used a simple sine function to model a mean temperature index throughout the year. 
(Temperature data from. 23 ) The mean air temperature index given by 

T a i r (t) = sin ( — — where t is in weeks from 17-Oct.-2010. 

w V 365.25 J ' 

This is then used to create a delay functions that varies from a minimum of r p j during 
summer and r p jy during winter 



T, 



V 



{ijpJii + TpJo) - (T pM - T pJo ) T air (t)) /2. 



In other words, infections result faster in the summer than in the winter months. This 
would be expected since warmer temperatures mean faster growth rates for bacteria and their 
invertebrate hosts as well as triggering emergence of sporesE32H 

In addition, there is a direct influence of temperature on the infection rate. This heat 
index is given by 

H(t) = l + kT air (t) . 
Finally, the infection rate is given by 



(3 (t) = u(t) a p H(t) P(t - t p , 6 } 



V) 
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2.6.2. Ouest. For the Ouest region we use the same formulation as in Artibonite, however, we 
found that tidal range appeared to significantly affect infection rates as well. The maximum 
tidal range each day (there are two) m (t)P^ is averaged over an interval 9 m so the running 
average precipitation rate is 

M(t - r m , m ) = — ^2m(t - r m - j) - M , 

m j=0 

where r m is the delay in the tide's affect on infection rate. Here, r m is fixed, the effect of 
temperature (water or air) on lengthening or shortening the response in infection rate was not 
found to be sufficient to warrant adding another function and additional parameters. Thus, 
the overall infection rate for Ouest is 

(t) = u(t) [a p H(t) Pit - Tp , 9 P ) + a m M(t - r m , 9 m )\ . 



3. RESULTS 

3.1. Parameter Fitting and Model Selection. The following Table [2] displays parameter 
values for each region obtained through curve-fitting to cumulative reported cases. 



parameter 


Artibonitc 


Ouest 


units 


description 


6 p 


4.61 (0.319) 


1.96 (0.526) 


weeks 


averaging window for precip. 


1~p_lo 


0.4019 (0.3261) 


3.800 (0.3757) 


weeks 


minimum delay for precip. effects 




7.1805 (0.2056) 


6.500 (0.4240) 


weeks 


maximum delay for precip. effects 


k 


0.5198 (0.0177) 


0.5812 (0.03050) 


o C -l 


temp. -precip. interaction level 


a p 


6.8806 (0.1690) x lO" 3 


1.1502 (0.08582) x lO" 3 


(mm Xweek) -1 


infection rate per mm rain 


dm 




1.00 (0.91) 


weeks 


averaging window for tidal range 


I'm 




1.960 (0.494) 


weeks 


delay for tidal range effects 


OL m 




2.8128 (1.3360) x 10" 4 


(cm Xweek) -1 


infection rate per cm tide range 


M 




24.791 (24.094) 


cm 


baseline of tidal range effect 


r 


0.06100 (0.00147) 


0.02610 (0.00139) 


week -1 


decrease in susceptibles per week 



TABLE 2. Parameter values obtained from model calibration. Numbers in paren- 
theses are 95% confidence intervals (CI's). 



After initially restricting the parameter search space by visually fitting the New Cases 
predicted to the New Cases data, we used the Berkeley-Madonna curve-fitting routine on 
cumulative cases and in an iterative fashion. We first optimized the averaging windows 
(0's) and delays (r's), then temperature interaction level (k), infection rates (en's), and the 
improvement rate (r). We then re-optimized the averaging windows and delays and then a 
final fit for the parameters a p , r and for Ouest a m . This was done for two reasons, first, the 
high variability of the data and the non-linearities of the averaging windows and delays lead 
to a number of local minima of the RMSD's, and second we wanted to search in a particular 
basin that corresponded to our visual fit of the New Cases. Allowing the software to search 
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the entire parameter space all at once would inevitability cause it to jump from basin to basin 
to find a global minimum for just cumulative cases. 

The Artibonite model with tide was not included in Table [2] since inclusion of tide did not 
improve the model (see Tables [3] and [1]). The statistics for the model fit are given in the 
following two tables. Table [3] is for cumulative cases predicted by the model compared to 
cumulative case data. 



All cases 


statistic 


Artibonite (tide) 


Artibonite (no tide) 


Ouest (tide) 


Ouest (no tide) 


data points 


77 


77 


74 


74 


parameters 


10 


6 


10 


6 


adj RMSD 


2181.40 


2121.78 


3708.59 


6336.53 


adj R 2 


0.9930 


0.9933 


0.9968 


0.9908 


AIC 


1193.20 


1185.40 


1271.40 


1352.56 



TABLE 3. Model predictions versus data statistics for cumulative number of cases: 
Root mean squared deviations (RMSD), Coefficient of determination (R 2 ), and 
Akaike information criterion (AIC). Degrees of freedom for the statistics are adjusted 
by the number of parameters fit in the calibration process. 



For the full model in either region (model including tides) the parameters m ,r m , M , and 
a m are added and the model is re-optimized. For Artibonite an F-test for the nested models 
gives the following results F = 0.043 d.f. = (4, 67) and the p value is 0.9964, indicating the 
inclusion of the extra detail had almost no effect and would be unjustified. For the Ouest 
region an F-test for the nested models gives the following results F = 33.63 d.f. = (4, 64) 
and the p value is 4.208 x 10~ 15 , indicating the tidal data significantly improved the model 
fit. 



New cases 


statistic 


Artibonite (tide) 


Artibonite (no tide) 


Ouest (tide) 


Ouest (no tide) 


data points 


77 


77 


74 


74 


parameters 


10 


6 


10 


6 


adj RMSD 


9.601 


9.348 


21.30 


23.08 


adj R 2 


0.7913 


0.8027 


0.3932 


0.3681 


AIC 


357.62 


349.97 


461.92 


470.28 



TABLE 4. Model predictions versus data statistics for new cases: Root mean squared 
deviations (RMSD), Coefficient of determination (i? 2 ), and Akaike information cri- 
terion (AIC). Degrees of freedom for the statistics are adjusted by the number of 
parameters fit in the calibration process. 
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Table [4] is for new cases predicted by the model compared to new case data. The statistics 
in Table [4] are done on the square roots of the values in order to produce more uniformly 
distributed data and stabilize the variance somewhat. The F-test for the Artibonite nested 
models using new cases gives the following results F = 0.0755 d.f. = (4, 67) and the p value 
is 0.9894. For Ouest the difference is again highly significant with F = 3.9599 d.f. = (4,64) 
and the p value is 0.0062. 

3.2. Lag times. The total delays in response to precipitation and tides are the sum of the 
averaging window and the delay function. For precipitation the minimum and maximum 
delays for Artibonite are 5.0 and 11.8 weeks, and for Ouest they are 5.7 and 8.5 weeks, 
respectively. The shorter delays are very similar in the two regions during the warmer months 
but during the cooler months the response time in Artibonite is 3 weeks longer. These long 
delays are similar in magnitude to delays reported from a study of Cholera in Zanzibar, East 
Africa (8 weeks fixed delay) and the shorter delays (4 weeks) to those in Bangladesh.^ For 
Ouest estimated delay from response to changes from tidal range was about 3 weeks. But 
since influence of tidal range has not been quantitatively reported elsewhere in the literature, 
we have nothing to compare this number to. 

Although, rainfall data was available only to week of 26-Feb-2012 (week 71), we run the sim- 
ulations for another month, until l-Apr-2012 (1st day of week 76), just short of the minimum 
time lag between precipitation events and new cases. 

3.3. Vaccination. A program to vaccinate the most at risk populations began in the 2nd 
week of April and ended in mid June. Each site (Ouest and Artibonite Dept) vaccinated 
about 50,000 persons, and each site had about 91% 2nd dose coverage. The administration of 
the first dose was staggered by age groups (beginning first with 10 year olds and up) because 
the Ministry of Health had a measles, rubella and polio vaccine catch-up campaign for children 
under 10 years of age that was taking place at the same time last April 2012^] 

In the Ouest Department, GHESKIO vaccinated adults, adolescents and children over 10 
years of age from April 12-23, 2012 and children under 10 from May 26 - June 3. Vaccine was 
provided to 99,877 persons (52,357 with the first dose, and 47,520 the second dose) living in 
the slums of Port-au-Prince and surrounding villages]^] 

In the Artibonite Department, PIH vaccinated 32183 people in rural Bocozel and 13185 
people in Grand Saline with 90.8% of those people confirmed to get the 2nd dose (or 41194 
for both locations). The campaign started April 15th 2012 and ran until June 10 2012. Here 
too, children under 9 year old were vaccinated in the second half of the time period because 
of the MMR and Polio vaccination campaignjjj] 

With these basic facts we constructed a crude vaccination schedule (Table [5J using the 
following assumptions: 

^Communicated by Jordan Tappero, MD, MPH (CDC/CGH/DGDDER, Atlanta, GA). Nov. 29, 2012. 
Communicated by Jean W. Pape, MD (GHESKIO, Weill Cornell Medical College, Port-au-Prince, Haiti). 
Nov. 29, 2012. 

"Communicated by Louise Ivers, MD (Partners in Health/ZL, Cange, Haiti). Nov. 29, 2012. 
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1) approximately 25% of the population is under 10 years old; 

2) the second dose was administered 14 days after the first dose was given^ 

3) the immune response took hold about 8.5 days after the second dose was given^ 

4) we used the average number of people vaccinated per day over a 12 day period for adults 
and 9 days for children. 



Ouest 


Adult I st dose 


39,268 


April 10 - April 23 


Adult 2 st dose 


35,640 


April 26 - May 7 


Adult immune response 


65 - 85% 


May 4 - May 15 


Child 1 st dose 


13,089 


May 26 - June 3 


Child 2 st dose 


11,880 


June 9 - June 17 


Child immune response 


65 - 85% 


June 17 - June 25 




Artibonite 


Adult I st dose 


34,026 


April 15 - April 26 


Adult 2 st dose 


30,896 


April 29 - May 10 


Adult immune response 


65 - 85% 


May 7 - May 18 


Child 1 st dose 


11,342 


May 19 - May 27 


Child 2 st dose 


10,298 


June 2 - June 10 


Child immune response 


65 - 85% 


June 10 - June 18 



TABLE 5. Simulated vaccination schedule for Artibonite and Ouest 



We ran simulations including vaccination by removing, in Artibonite, 2574.67 at risk persons 
over age 10 per day starting on 7-May-2012 and ending on 18-May-2012 (30896 total), then 
removing 1144.22 at risk children per day starting on 10-June-2012 and ending on 18-June- 
2012 (10298 total). In Ouest we removed, 2970 at risk persons over age 10 per day starting 
on 4-May-2012 and ending on 15-May-2012 (35640 total), then removing 1320 at risk children 
per day starting on 17-June-2012 and ending on 25-June-2012 (11880 total). 

Since the efficacy of the vaccine (oral Shanchol) is between 65 and 85%P^we did three runs 
65, 75, and85% for each department. The 75% run was our mean. To get the bounds we 
added the 95% prediction interval to the 65% run and subtracted it from the 85% run. 

3.4. Simulations and Projections. 

3.4.1. Predictions compared to observations for cumulative and new cases. The following Fig- 
ures [I] and [2] show the model predictions compared to observations for cumulative number 
of cases in Artibonite and Ouest regions, respectively. Similarly Figures [3] and [4] show the 
model predictions compared to observations for new number of cases in Artibonite and Ouest 
regions, respectively. Prediction intervals were calculated only for the cumulative numbers 
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since the final model fitting was done on these numbers. Confidence intervals are shown for 
incidence. The match for the trends in new cases match fairly well, the slope of the expected 
(model) regressed against observed (data) is nearly one in both departments (see Figures [5] 
and [6]) even though there is a substantial amount of unexplained variance. Whether this is 
due to the crude spatial resolution or other factors remains to be seen. 

Prediction intervals (Pi's) and confidence intervals (CPs) were calculated using the delta 
method adapted for differential equations (see, for example, Ramsay et al. 2 ^). After 30-Sept- 
2012 rainfall data from NASA was unavailable, for simulations after that date we did 13 runs 
using rainfall patterns from each of the 13 prior years. The mean of those simulations was 
used and the variance of the 13 runs, at each time step, was added to the variance from the 
estimation procedure before computing the Pi's. 
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Predicted and Reported Cumulative Number of Cases vs. Time 
Antibonite (17-Oct-10 to l-April-12 with projections to 27-Feb-13) 

140000 t 1 1 1 1 1 1 1 1 1 




Figure 1. Artibonite. The predicted cumulative number of symptomatic indi- 
viduals, against total reported cases to l-Apr-2012. Projections are from then 
to end of February. Projections using the vaccination schedule (red) begin on 7- 
May-2012, approximately three weeks after beginning the vaccination program 
in Artibonite. All projections after ll-Nov-2012 are based on runs using prior 
13 years of precipitation, and Pi's include the variance of those data (see text). 
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Predicted and Reported Cumulative Number of Cases vs Time 
Ouest (17-Oct-10 to 1 -April- 12 with projections to 27-Feb-13) 

300000 1 1 1 1 1 1 1 1 




Figure 2. Ouest. The predicted cumulative number of symptomatic individ- 
uals, against total reported cases to l-Apr-2012. Projections are from then to 
end of February. Projections using the vaccination schedule (red) begin on 4- 
May-2012, approximately three weeks after beginning the vaccination program 
in Ouest. All projections after ll-Nov-2012 are based on runs using prior 13 
years of precipitation, and Pi's include the variance of those data (see text). 
Note that for the Ouest region, the model begins at the fourth week. We assume 
that the low initial numbers in the first three weeks are a result of immigration 
of cases from the Artibonite region. The model therefore uses data for the first 
four weeks - assumed immigration numbers for the first three weeks and the 
initialization of the model from data for the fourth week. 
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Predicted and Reported Number of New Cases per Week vs Time 
Artibonite (17-Oct-10 to 27-Feb-13) 




17-Oct-10 12-Dec-10 6-Feb-ll 



1-Sep-ll 13-Nov-ll S-Jan-12 4-Mar-12 29-Apr-12 24-Jun-12 19-Aug-12 14-Oct-12 
Date 



9-Dec-12 3-Feb-13 



Figure 3. Artibonite. The new symptomatic individuals, vs. time. Circles - 
observed; solid line - model prediction; Dashed lines - 95 percentile confidence 
intervals for model projections. Red line is projections with 75 percent vaccine 
efficacy. 
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Predicted and Reported Number of New Cases per Week vs Time 
Ouest (17-Oct-10 to 27-Feb-13) 

16000 -i — — i-s 1 — — I — — i — — i — — i — — I — — i — — i — — i — — I — 




17-Oct-10 12-Dec-10 6-Feb-l 1 3-Apr-ll 29-May-ll 24-Jul-ll 18-Sep-ll 13-Nov-ll 8-Jan-12 4-Mar-12 29-Apr-12 24-Jun-12 19-Aug-12 14-Oct-12 9-Dec-12 3-Feb-13 

Date 



FIGURE 4. Ouest. The new symptomatic individuals, vs. time. Circles - observed; 
solid line - model prediction; Dashed lines - 5, 50, and 95 percentiles for model 
projections based on past 13 years precipitation records. 
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Predicted vs Reported number of New Cases per Week 
Artibonite (17-Oct-10 to l-April-12) 
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Figure 5. Artibonite. The predicted new symptomatic individuals, against 
weekly reported cases to l-Apr-2012 (square root transformed). A regression 
line matching the main diagonal (45°) dashed line would show an optimal fit, 
the discrepancy is due in part to fitting on the cumulative numbers. 
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Predicted vs Reported number of New Cases per Week 
Quest (7-Nov-lO to 1 -April- 12) 
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FIGURE 6. Ouest. The predicted new symptomatic individuals, against weekly re- 
ported cases to l-Apr-2012 (square root transformed). A regression line matching 
the 45° line would show an optimal fit, the discrepancy is due in part to fitting on 
the cumulative numbers. 
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3.4.2. Epidemic projections for Artibonite. The model projects that by the end of February 
2013, Artibonite would have seen between 117 and 130 thousand cholera cases without vaccine, 
and between 110 and 124 thousand with the implemented vaccination program (a decrease 
in about 7 thousand cases; see Figure [7]). The model also predicted that about 15002 people 
would have gotten ill between the onset of the immune response (6-May) and the end of the 
simulation (27-Feb) this represents about a 47% decrease in the number of people that would 
have gotten cholera. The F value comparing the mean values of two models to the actual cases 
between 6-May and 2-Dec is 9.88 with 25 degrees of freedom for both models, yields a p- value 
of 1.06 x 10~ 7 . This suggests that the different models (with and without vaccination) show a 
substantial difference in how well they match cumulative number of cases data. There seems 
to be better agreement between the projected and observed incidence as well (see Figure [8]) 
where the number of new cases is well below what would have occurred without vaccination. 
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Projected Cumulative Cases 
Antibonite (to week of 27-Feb-13) 
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FIGURE 7. Artibonite. The projected total symptomatic individuals, vs. time. 
Circles - observed; solid line - model prediction; Dashed lines - 5, 50, and 95 percentiles 
for model projections based on past 13 years precipitation records. Red lines are 
model runs with vaccination schedule that occurred in Spring 2012 (see text). 
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Predicted and Reported Number of New Cases per Week vs Time 
Artibonite (l-Jan-12 to 27-Feb-13) 




l-Jan-12 29-Jan-12 26-Feb-12 25-Mar-12 22-Apr-12 20-May-12 17-Jun-12 15-M-12 12-Aug-12 9-Sep-12 7-Oct-12 4-Nov-12 2-Dec-12 30-Dec-12 27-Jan-13 24-Feb-13 

Date 



FIGURE 8. Artibonite. The projected new symptomatic individuals, vs. time. Cir- 
cles - observed; solid line - model prediction; Dashed lines - 95 percentile confidence 
intervals for model projections. Red line is projections with 75 percent vaccine effi- 
cacy. 
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3.4.3. Epidemic projections for Ouest. For Ouest, the model projects that by the end of 
February 2013, Ouest would have seen between 244 and 266 thousand cholera cases without 
vaccine, and between 241 and 264 thousand with the implemented vaccination program (a 
decrease in about 3 thousand cases; see Figure 9). The model also predicted that about 33909 
people would have gotten ill between the onset of the immune response (6-May) and the end 
of the simulation (27-Feb) this represents about a 8.4% decrease in the number of people 
that would have gotten cholera. The F value comparing the mean values of two models 
to the actual cases between 6-May and 2-Dec is 1.47 with 21 degrees of freedom for both 
models, yields a p- value of 0.17. This suggests that the different models (with and without 
vaccination) do not substantially differ in how well they match cumulative number of cases 
data. Note that this does not mean that the vaccination program did not have an impact, 
but only that the model at this point in time could not distinguish between a scenario with 
and without vaccination. There was a slowdown in number of new cases in starting in mid 
February and running through March which dropped the cumulative case number below the 
mean model projections (but still within the prediction interval), however, this occurred before 
the vaccination program took place. There was another smaller slowdown in July but this 
was being made up for by a rapid growth in numbers in November (see Figure [9]). It would 
seem that these oscillations are in spite of, rather than because of vaccination, since the model 
predicts a steady, albeit slow, additional decline in new cases(see Figure 10). 
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Projected Cumulative Cases 
Ouest (to week of 27-Feb- 13) 

270,000 t 1 1 1 1 




FIGURE 9. Ouest. The projected total symptomatic individuals, vs. time. Circles 
- observed; solid line - model prediction; Dashed lines - 5, 50, and 95 percentiles for 
model projections based on past 13 years precipitation records. Red lines are model 
runs with vaccination schedule that occurred in Spring 2012 (see text). The slow 
down in growth in late February and March was before the vaccination program 
began. 
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FIGURE 10. Ouest. The projected new symptomatic individuals, vs. time. Circles 
- observed; solid line - model prediction; Dashed lines - 5, 50, and 95 percentiles for 
model projections based on past 13 years precipitation records. 
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4. Discussion 

Use of environmental data to model the dynamics of cholera in Haiti has been hampered 
by the lack of historical records. We use NASA satellite data to address this problem. This 
study shows that with environmental data of sufficient detail and quality, projections of disease 
progression can be made with sufficient lead time to prepare for outbreaks. The lag times 
of over five weeks means that if even rudimentary but reliable meteorological and coastal 
records are kept, preparations and resources can be more focused. The gathering of basic 
weather information is simple and inexpensive and should be made standard procedure when 
any agency takes part in interventions, particularly when the environmental component of 
the epidemiology is so well established. 

In addition we explored the hypothesis that, at least in the Ouest region, tidal influences 
play a significant role in the dynamics of the disease. It appeared that tidal range rather than 
the height of the tide itself had the strongest influence. Some connection to tidal influences 
should be expected where large populations are in close contact with bays and estuaries, and 
humans are consuming local seafood JHHH1 

We also affirmed the longer time lags (8 - 10 weeks) found in previous studies from Africa^ 
and shorter ones (4 weeks) in Bangladesh.^ 

Over the course of the epidemic the incidence has been tapering off. There has been steady 
and continued effort to improve hygiene and living conditions, however, the areas where the 
greatest strides are made are those where people leave the camps to return to normal living 
conditions and employment. The declining numbers of those at risk in the overall population 
belie the fact that many local populations are still without basic hygienic facilities. This was 
reflected in the model by setting a level uq to 5% of the original number of susceptible. This 
value matches the approximate 5% of the population that still remain displaced after the 2010 
earthquake.^ 

On top of predicting when and how many cholera cases will rise with Haiti's climate and 
tides, any modeling to predict the effectiveness of interventions (such as vaccination) should 
consider these patterns. Considering that cholera may be maintained in the environment 
outside the human chain of infection is essential to planing effective prophylaxes and inter- 
ventions. 

Using these models we were able to assess, to some degree, the relative effectiveness of the 
recent vaccination program in Artibonite and Ouest. The discrepancy between the apparent 
effectiveness of vaccination in the two regions is perhaps not that puzzling when one considers 
the number vaccinated relative to the size of at risk population. In Artibonite about 41 
thousand people and in Ouest about 47 thousand people received both doses of the vaccine. 
However, our model suggests that in Artibonite the at risk population by 6-May was about 
61 thousand whereas in Ouest it was still over 428 thousand, seven times the number in 
Artibonite. Thus in Artibonite 67% of the most at risk population apparently received the 
vaccine, while in Ouest only 11% did. Further, in both regions 100 people receiving vaccination 
does not mean 100 people protected. The vaccine is about 75% effective, and only 21% of 



26 



THE BIOCLIMATE OF CHOLERA: THE CASE OF HAITI 



people who are infected with cholera show symptoms, 0.75 x 0.21 = 0.1575. So as a rough 
calculation we might expect 0.1575 x 41,000 = 6457.5 and 0.1575 x 47,000 = 7402.5 people 
protected directly and additional due to the reduction in environmental loading of cholera 
bacteria and the force of infection. Ultimately, more people will be protected in Ouest but it 
will take longer to become apparent due to this dilution effect. 

4.1. Conclusion. Although progress has been made in the past two years in modeling the 
Haitian cholera epidemic, only one model, so far, has accounted for seasonal variation due to 
precipitation, and none have examined temperature or tidal influences. Our paper shows that 
with basic climatic and tidal records and a relatively unsophisticated modeling procedure, 
not only can these factors be accounted for but also the time lags between climatic events 
and outbreaks can be identified. This approach looks at specific climatic events on the scale 
of a week rather than just seasonal patterns. We use daily tidal range as a predictive factor 
for cholera epidemics for the first time in a modeling paper, and we use real time daily pre- 
cipitation estimates. These provide a level of detail of environmental events in conjunction 
with the lag times and estimates of core population size that can help evaluate intervention 
(such as vaccination) and public hygiene efforts. In order to show this we examined recent 
vaccination efforts in the Haitian cholera epidemic. Complex environmental patterns incor- 
porated in epidemic models allow us remove a large source of variability and bring into relief 
intervention efforts by identifying deviations from the current. 
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